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In this paper, I study the oxygen isotope effect (OIE) in cuprates. I introduce a simple model that 
can explain experiments both qualitatively and quantitatively. In this theory, isotope substitution 
only affects the superfiuid density, but not the pseudo-gap. Within the spin-charge separation 
picture, I argue that the spinon-phonon interaction is in the adiabatic limit, and therefore within the 
Migdal-Eliashberg theory, there is no isotope effect in the spinon mass renormalization. On the other 
hand, I show that the holon-phonon interaction is in the non-adiabatic limit. Therefore, the small 
polaron picture is applicable and there is a large mass enhancement in an isotope-dependent way. 
Our theory explains why upon 16 0/ 1S substitution, the superconducting transition temperature 
changes only in underdoped cuprates, while there is no considerable OIE at the optimal doped 
as well as the overdoped cuprates. Additionally, in contrast to the conventional superconductors, 
we obtain OIE on the superfiuid density for whole superconducting region in agreement with the 
experimental observations. 



I. INTRODUCTION 

Finding the underlying mechanism of the high tem- 
perature superconductivity in cuprates [l| is one of the 
most challenging and outstanding problems in theoreti- 
cal physics. Observation of strong isotope effect on both 
the transition temperature and the superfiuid density in 
cuprates , indicates the importance of the electron- 
phonon interaction in understanding the physics of high 
Tc [7|, |8| . Experiments show strong OIE on T c only in 
underdoped cuprates (see Fig. 1), while there is no con- 
siderable OIE in overdoped cuprates. On the other hand 
OIE on the London penetration depth (A a & (0)) (in-plane 
penetration depth) has been reported for both sides 0- 
fl2j . In both cases the isotope exponent decreases as 
we approach the optimal doping from the underdoped 
side. There is also an unusual correlation between iso- 
tope effect on T c and A a6 (0) [3 (See Fig. 2). It is im- 
possible to explain such effects using BCS theory [l4| 
or its extensions such as Migdal-Eliashberg [ll| theory. 
The reason is that within BCS theory which is based 
on the adiabatic electron-phonon approximation (ij, the 
isotope exponent of T c which is defined as a — — ^^j, 
is around 1/2 14]. On the other hand, the electron ef- 
fective mass is m* — m (1 + A) [16j |, where A is the di- 
mensionless phonon mediated attraction coupling, and is 
isotope independent. Therefore, theories based on the 
adiabatic approximation, predict absence of the isotope 
effect on the superfiuid density (n s /m*) which is incon- 
sistent with experiments. Therefore, we have to look for 
other theories of superconductors. In this paper we start 
from Anderson theory of high Tc [|| which is based on 
the strong correlation physics and the spin fluctuation 
pairing mechanism. This theory is very successful in ex- 
plaining many aspects of cuprates but it does not take 
phonons into account. It is very important to explain the 
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oxygen isotope effect within this successful theory. Here 
we use this general framework and then by adding elec- 
tron phonon interaction, it is shown that the observed 
unusual isotope effects on T c as well as A a f, (0) can be 
explained. 

We start from t-J-Holstcin model as our model Hamil- 
tonian. In our model for simplicity, electrons are only 
coupled to a single Einstein phonon mode (w E ). Then 
we use the spin-charge separation picture and we imple- 
ment it by the slave boson method. Our idea is as fol- 
lows: in underdoped cuprates, the superconducting tran- 
sition temperature is controlled by the superfiuid density 
(n s /m*), which is mostly determined by that of holons, in 
our theory. On the other hand in overdoped cuprates, T c 
is controlled by the pseudo-gap which is equal to pairing 
order parameter of spinons A s in our theory. Electron- 
phonon interaction affects the superfiuid by holon mass 
renormalization in an isotope-dependent way, while it 
does not affect the pseudo-gap much. In our theory, 
we argue that holons are strongly coupled to phonons 
and therefore, we are in the non-adiabatic limit, while 
spinons are weakly coupled to phonons and we should use 
the adiabatic limit calculation (Migdal-Eliashberg the- 
ory) for them. Therefore, < uj e < E s , where E\, 
and E s are typical energies of holons and spinons respec- 
tively. It is also shown that in the presence of spinon pair- 
ing (pseudo-gap), the holon-phonon coupling constant^, 
renormalizes to A s 7. These facts together can explain 
the observed OIE on T c as well as on A -2 (0) cx ^f-. 

II. METHOD 

Let us start from the t-j[l7i[II| model which is defined 
as: 

Ht-j = -* £ PgcI^Pg + JJ2^-S, (1) 

(i,j),ir id 

where Pq is the Gutzwiller projection operator which 
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FIG. 1: Measured oxygen isotope coefficient of T c (a) for 
doped YBCO as a function of the reduced transition temper- 
ature (T c = T c /T Cl max) from different samples. After[3] 




FIG. 2: A color plot of the OIE shift AX ab (0) /Kb (0) 
versus the OIE shift ~AT C /T C for Yi_ 2 ,Pr a; Ba2Cu307-i, 
YBa2Cu408, and La2-^SrxCu04. Squares are the fiSK data 
obtained in the present study. Circles are bulk /iSR data for 
Yi_ a; Pr 2 ,Ba2Cu307-5 (Refs. 4 and 12 ) and LE /xSR data for 
optimally doped Yi_ a; Pr :E Ba2Cu307_,5 (Ref. 10 ). Triangles 
are torque magnetization and Meissner fraction data for La2- 
xSrxCu04 (Refs. 6 and 9). Different lines are discussed in 
Ref. 13 



- ]T 2 J A a (4/J 4 - f}jl) + h.c. (3) 

<i,i> 

where the following notations have been used: 

Mm) = ^(44 "44) ( 6 ) 

Form the t-J model one may expect that J = J/4 but 
in literature it has been discussed that the best choice 
is J = 3J/8[I3|- This model has been studied by many 
authors in the past three decades. It is well known that 
this model can lead to the d-wave pairing of spinons pjj, 
i.e. A s (±x) = A s and A s (±y) = — A s . 

Now let us consider the electron-phonon interaction. 
To do so we add the following Hamiltonian to the t-J 
model. 

H' = ^uJ E a\a q + ^2 7fc, 9 (a q + al q ^ c\ +q(J c k ^ (7) 

q k,q 

where Ck.a is the Fourier components of PgCi^Pg = 
fi.ahj. If we treat fcrmionic operator, then since 

the typical energy of electrons is around J and is much 
larger than u E , so that we can apply the BCS theory, 
or the Migdal-Eliashberg theory as its extension. BCS 
approximation lead to the following pairing term: 

- Vk , k ' < A'.^-k'.x > e -fc4Cfc,t ( 8 ) 

k,k' 

Now if we translate everything to the slave boson lan- 
guage in real space, within mean-field approximation we 
can substitute Ci^Cj^ = fi,lfj,fh\hj by the following 
form: 



removes doubly occupied states. Within the slave bo- 
son formalism, electrons can be decomposed as cj a — 
fi a hi along with the physical constraint on each site: 

h\hi + J2a ficrfi,cr = 1 which implements the Gutzwiller 
projection. Now if we decouple spinons (spin sector) from 
holons (charge sector), we can write two following effec- 
tive Hamiltonians for each sector: 

H h = - J2 t**^ ~ ^ h t hi ( 2 ) 

<i,j> i 



A ft (i,j)fi,if M + A s (i,j)h]h] - A s (i,j)A h (9) 

where Ah = (hihj). Now assuming a very short 
range interaction, we obtain the following effective inter- 
action: 

H>_ s = -VAl ]T A s 4 t /j 4 + h.c. (10) 

<i,j> 

H' h _ h = -VAl ]T ^h\h) + h.c. (11) 

<i,j> 

where A/, ~ p (doping). Let us assume that VA^ ~ 
Vp 2 <C J, so we can neglect this phonon mediated pair- 
ing term and therefore, the d-Wave nature of the spinons 
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does not change. From the above we see that the cou- 
pling constant of phonon mediated spinon spinon attrac- 
tion is renormalized by A 2 factor and that of holons by 
A 2 due to strong correlation effects. It is easy to show 
that, V oc 7 2 , therefore, we can interpret these renor- 
malization factors as the renormalization of the coupling 
constant of spinon-phonon interaction to Ahjk , q and that 
of holon-phonon interaction to Ah"fk, q - Therefore, we can 
substitute the electron phonon interaction term by sum 
the following two terms: 



H s ^ ph = &hlk, q (a, + al q ) fl +qt J k ,„ (12) 

k,q 

H h -p h = &slk, q (a q + al q J h\ +q h k (13) 

k, q 

Now let us consider spinon-phonon interaction. The 
typical energy of spinons is around J. J is around 1500 
Kelvin while the typical energy of optical phonons (lu e ) 
is a few hundred kelvin. Therefore, spinon-phonon in- 
teraction is in the adiabatic limit and therefore, we can 
apply Migdal-Eliashberg theory. In this regime the mass 
renormalizes as: m* = m (1 + A), where: 



— a 2 (w)F(w) 



(14) 



From our simple model it can be shown that A oc jj^r 



and since uj„ oc 



= , A is isotope independent. There- 



fore, oxygen isotope substitution does not enhance the 
effective mass of spinons. This agrees with experiment 
where there is no isotope effect on Fermi velo city by Laser 
ARPES, while there is shift in kink energy [19j . 

On the other hand, holons are hard-core bosons. They 
usually condense at the bottom of their energy band. 
So their effective band- width is much smaller than their 
actual band- width. To have a better idea, for the moment 
let us assume that they are fermions. In that case their 
Fermi energy will be around 4x ftp , where p is doping. 
Therefore, their effective bandwidth is at most of the 
order of tpx / ■ Now if w E is larger than the typical kinetic 
energy of holons, the electron phonon interaction is in the 
non-adiabatic limit. In summary if tpxf < uj e < J then 
spinon-phonon interaction is in the adiabatic limit and 
holon-phonon interaction is in the non-adiabatic limit. 
For this limit it is easier to rewrite Hh + H v h + Hh-ph in 
the following way: 



co E a\a q 



<i,j> 



\— ^ s in, q h{h n (di q +d q yi5) 




FIG. 3: The oxygen isotope effect on T c (a) vs. T c /T c ^ max . 
Red square are the data extracted from Fig. 1 and the blue 
line is our theoretical calculation 



a single polaron theory [20, Hl(. Therefore, we 
can do the powerful Lang-Firsov transformation [22j], 
i.e. H -> H = e~ s He s where: S = 
* Ey /2NMu„ ^, q A s hih n (dl q - d q ) . As an ap- 

oxi 

Therefore, we finally have: 



exp(ik.R n ) 

proximation let us assume that < j n<q >— jq ^= — -. 



H ~ {Atx s (cos k x + cos k y ) + fi h ) h\h k 

k 

+5> £ (4^ + 1/2) 

9 

hq 3 h— qi —q 2 —q 3 (16) 

91.92,93 

where we have used the following notations: 



Xs (T) 



= P -g 2 (T) 



Xs (T) 



7 2 (T) ^7oA?(T) 



4tAn = V 



MujI 



(17) 
(18) 

(19) 



in which M is the ion mass. Now if we expand the 
energy of holons around k = we have: 



e h (A:) = 
m *h (T) 



k 2 



2mt 



Vh 



2f s (T) 



e 9 {T) m h 



(20) 
(21) 



q ,n.a 



Oxygen isotope effect on T c 



The above model has been extensively studied. In 
the non-adiabatic regime, we can use the results for 



The oxygen isotope effect on T c is determined by the 
a isotope exponent which is defined as: 
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Q 



M dT c lu. dT c 



o 



2 T r dui c 



(22) 



In the underdoped region, T c is determined by the 
Bose-Einstein condensation (BEC) transition tempera- 
ture of holons. If we assume holons as free 2D bose gas, 
then the BEC transition temperature is: 



T c — Tbec 



2-Ktp 
m* h (T c ) 



(23) 



Therefore, we have: 



= 1 u E dm* h (T c ) 
2m* h (T c ) <kj E 

Since m* h — e B °mh and g oc we have: 



a = ( T c) 
From the definition of g 2 , we have: 



(24) 



(25) 



a 



2tX A 2 s (t c ; 



UJ, 



(26) 



Since A s (T c ) is a decreasing function of doping and 
at the optimal doping it becomes very small, isotope ex- 
ponent dies off as we approach the optimal doping and 
it finally becomes negligible. On the overdoped cuprates 
however, pseudo-gap controls T c and as we discussed be- 
fore, if Vp 2 — 4t\op 2 <C J, electron phonon interaction 
does not affect pseudo-gap much and therefore, we do not 
expect isotope effect in overdoped side of the supercon- 
ducting phase diagram. In our numerical calculation we 
have studied the phase diagram of the t-J for J = t/3. 
We have also chosen Ao = 3u E /t ratio. We obtain a very 
good fitting between our theoretical calculations and the 
experimental data (see Fig. 3). 



Oxygen isotope effect on the London penetration 
depth 

The oxygen isotope effect on the superfluid density or 
equivalently on the London penetration depth, is defined 
as: 



a _ Mp dX ab _ 1 



X ab dM Q 



u E dXj 



4X- 2 du E 



(27) 



According to the Ioffe-Larkin formula [23| the physical 
superfluid density is related to the superfluid density of 
spinons and holons in the following way: 




FIG. 4: The OIE on T c , a (dashed blue line) and the OIE on 
the London penetration depth, |/3| (black line) versus doping. 
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FIG. 5: The OIE shift A ,Vffl versus the OIE 

A ab(0,) 



(P er - 

cent). Red squares are extracted data from Fig. 2. The black 
line is our theoretical calculation and the dashed blue line is 
the shifted black line. OIE on the thickness of the supercon- 
ducting sheet can cause this shifting. 



-l -l , -l 
Pph = Ph + Pb 



(28) 



Since condensation fraction of holons and spinons at 
zero temperature are p and 1 — p respectively, and from 
p = we have: 

i 771 



p^ (o) = ^ 



(29) 



p 1-p 

For small values of p, we have p p h (0) ~ -Ar. On the 
other hand X~ 2 (0) = ^§-p p h- Therefore, we have: 



_ 1 uj e dm* h 
4 to* du E 



5 2 (0) 



(30) 



Therefore, we have: 
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/3(r = 0) = _!AoA|(0) 

For the whole superconducting region, A s (T = 0) > 0, 
therefore, |/3| is always nonzero though it is a decreasing 
function of the doping. Despite the fact that both a and 
P depend on A 2 and we may expect |/3| = 0.5a, but 
we should note that they depend on A s at two different 
temperatures (see Figs. 4 and 5). From Eqs. (26) and 

(31), we have: \P\ /a = 0.5 (37$^) > 0.5. 

Discussion. — Our theoretical curve for A x A "^ffl ver- 

sus — has two important features. Firstly it has a 
nonzero inception since at the optimal doping, T c does 
not change upon isotope substitution, but A a fc (0) does. 
Secondly it has the same slope(0.42) as the empirical data 
but our curve is shifted down by a constant amount (see 
Fig. 5). In literature the cause of this shift has been 
discussed [T3|. One possible scenario is the change in 
the thickness of the superconducting sheet d Sl due to 
16q/18q substitution. Note that the 2D density of holons 



is related to the 3D one by: n 2D — n. lD d s . Thus A ab 2 oc 
and finally we have: = 0.5 - . 

If we assume 0.5-^^ = —1.6 then our theoretical curve 

d s 

fits the experimental data very well. So our theory pre- 
dicts OIE shift of the lattice spacing in the z direction. 

Conclusion. — I have studied the oxygen isotope ef- 
fect in cuprates. I have shown that within the spin-charge 
separation formalism, the unusual isotope effect on T c as 
well as the superfluid density can be explained by t-J- 
Holstein model. It is shown that 16 0/ 18 substitution 
only affects the superfluid density and it does not affect 
the value of the pseudo-gap energy. This model also pre- 
dicts the OIE on the thickness of the superconducting 
sheet. 
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